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Orbital dynamics: The origin of the anomalous optical spectra in ferromagnetic 

manganites 

R. Kilian and G. KhaHullinS 
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(June 30, 1998) 

We discuss the role ol orbital degeneracy in the transport properties of perovskite manganites, 
focusing in particular on the optical conductivity in the metallic ferromagnetic phase at low tem- 
peratures. Orbital degeneracy and strong correlations are described by an orbital t-J model which 
we treat in a slave-boson approach. Employing the memory- function formalism we calculate the 
optical conductivity, which is found to exhibit a broad incoherent component extending up to bare 
bandwidth accompanied by a strong suppression of the Drude weight. Further, we calculate the 
constant of T-linear specific heat. Our results are in overall agreement with experiment and suggest 
low-energy orbital fluctuations as the origin of the strongly correlated nature of the metallic phase 
of manganites. 



PACS number(s): 72.80.Ga, 71.27. +a, 71.30. +h, 

Perovskite manganese oxides, Ri^ x A x MnOs (where R 
and A represent rare-earth and divalent metal ions, re- 
spectively) exhibit rich physical behaviorEJ whose origin 
lies in the mutual interplay between spin, charge, orbital, 
and lattice degrees of freedom as well as in the strong 
correlations among electrons. The key elements of elec- 
tronic structure are the Mn 3d orbitals, which due to 
crystal-field splitting and strong Hund coupling are in 
a (< : 2g) 3 (eg) 1 configuration for x — 0, with holes being 
induced in the twofold degenerate e g orbitals for x > 0. 
Electrons in the ti g and e g orbitals are usually considered 
as localized and itinerant, respectively. 

Early work on the transport properties of manganites 
has focused primarily on the double-exchange scenario ,□ 
which also, competing with antiferromagnetism, results 
in a complex magnetic phase diagram. The renewed in- 
terest in manganese oxides sparked by the recent discov- 
ery of a very large magnetoresistance has led to a contro- 
versial discussion of the role of lattice and orbital degrees 
of freedom; the effect of onsite correlations among e g elec- 
trons, however, has mostly been discarded. On the ex- 
perimental side, these correlations seem to be strong even 
in the metallic state of manganites, despite the fact that 
double exchange strongly suppresses the spin dynamics 
below T c . This can be deduced, e.g., from the obser- 
vation of only a small discontinuity at the Fermi level 
in photoemission spectraJ3u Particularly striking in this 
respect are recent measurements of the optical conductiv- 
ity of Lai_ x Sr s Mn03.El The conductivity in the metallic 
state was found to consist of two components: (1) a nar- 
row Drude peak at lu < 0.02 eV with strongly suppressed 
weight and (2) a broad incoherent part extending up to 
w rj 1 cV. It is remarkable that the incoherent com- 
ponent remains finite even at T < T c where ferromag- 
netic moments have completely saturated, clearly indi- 
cating the presence of other degrees of freedom besides 
spin fluctuations. The fact that the incoherent part of 
the spectrum extends to rather high frequencies suggests 
these degrees of freedom to be of electronic origin, pos- 
sibly resulting from the orbital degeneracy of e g orbitals 
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as has been proposed by severaL.authors.il] 

Shiba, Shiina, and TakahashiQ ascribe the excitations 
leading to the incoherent structure of the optical con- 
ductivity to transitions between two fermionic bands 
within the degenerate e g orbitals, which they find to 
result in a spectrum extending up to fermionic band- 
width. While their model indicates the role of orbital 
degeneracy, it neglects the e g -electron correlations. A 
more elaborate treatment of both, orbital degeneracy and 
electron correlations, was done by Ishihara, Yamanaka, 
and Nagaosa.B Starting from an orbital-charge separation 
scheme these authors showed that orbital fluctuations are 
strong enough to prevent orbital ordering from develop- 
ing. Performing numerical calculations based on a phe- 
nomenological model that simulates orbital disorder by 
static randomness in the hopping matrix elements they 
further were able to recover the characteristic features 
of the optical conductivity spectrum, though no Drude 
component was obtained. 

In this paper we report on the first microscopic theory 
of the optical conductivity which combines strong cor- 
relations and orbital degeneracy. The transport prop- 
erties of manganites are shown to be highly incoher- 
ent even in the ferromagnetic phase due to strong scat- 
tering of charge carriers on dynamical orbital fluctua- 
tions. This gives rise to a broad optical absorption spec- 
trum extending up to bare bandwidth accompanied by a 
strong suppression of the Drude weight. The theory fur- 
ther accounts far the small values of specific heat found 
experimentally.!] 

Assuming complete ferromagnetic saturation of elec- 
tronic spins, which will henceforth be neglected, our 
starting point is an orbital t-J model: 



(1) 



(ij)-y a P 



where (ij) 7 indicates summation over manganese 
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nearest-neighbor bonds in spatial direction 7 S {x, y, z}. 
The operator c\ a = c\ a (1 — creates an electron at an 
empty site i with orbital pscudospin a for which we us 
the notation |= d 3z 2_ r 2 and 1= d x 2_-x. The anisotropic 
transfer matrix elements are given bylij 
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a/3 
x/y 
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allowing for interorbit transitions in the xy plane. The 
strength of the pseudospin interaction is controlled by the 
superexchange coupling constant J = 2t 2 /U, where U 
the onsite repulsion between spin-parallel e g electrons. 
Since U is believed to be large we assume J < t. The 
pseudospin operators of Eq. (|lj) are defined as 

~-X I y 1 ( z 1 /o _:r\ ^z ^ z 

4 V / ' 2 

with Pauli matrices cr 31 and <r z . We note that in the 
real system Jahn- Teller coupling further contributes to 
the J term of Hamiltonian ([I]) . Three-site hopping_terms 
that are contained in a general orbital t-J modefij are 
neglected for simplicity. 

We calculate the optical conductivity <t(uj) using the 
memory-function formalism.E3 While not rigorous this 
method yields exactly the leading terms of a high- 
frequency expansion of a(u) and is believed to give rea- 
sonably accurate results over the whole frequency range 
if no critical low-energy modes as in one-dimensional sys- 
tems exist. Within this framework the optical conductiv- 
ity <j(lo) is expressed via the memory function M(us) 



= Xo 



u + M (w) ' 



(2) 



where xo is the zero-frequency current-current correla- 
tion function. The memory function is given by 



M(w) 



1 

^Xo 



[/M-/(o)], 



(3) 



with f(uj) being the correlation function of force opera- 
tors F z — [J z , H] . The transport properties of the system 
are isotropic in the orbital-disordered phase, allowing us 
to select the current along the z direction which is of the 
simple form 



Jz 



iet 
~Nq3 



£(4t 



H.c. 
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where N is the number of lattice sites and a denotes the 
lattice constant. The force operator consists of two parts 
F z and Fj corresponding to the t and J term of Hamil- 
tonian ([j]) . Due to the relative smallness of the coupling 
constant J as compared to t the latter term is expected 
to contribute only little to the optical conductivity and 
will therefore be neglected. For the former term one finds 
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x {B ia cl +S ^c i+Jt p + H.c.) , 
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with the bosonic operator 



{2-n % + erf) <S qT + (erf - ia\ ) 5 al , 



where double occupancy of sites is excluded. 

We employ a slave-boson theory to express Hamilto- 
nian (Q) and the force operator, Eq. (^), in terms of 
fermionic spinon operators ft a and bosonic holon oper- 
ators hi. The pseudoparticles fulfill the local constraint 
/if /i| + fiifil + hi hi = 1 which is relaxed to a global 
one. On a mean-field level, the dynamics of spinous 
and holons are given by decoupling Hamiltonian (|l|) into 
fermionic and bosonic parts, introducing mean-field pa- 
rameters x — {fl-\fj])z and \fx — (hi), where x is the 
concentration of doped holes. The mean-field Hamilto- 
nian describes two fermionic Gutzwiller bands of width 
W s — 6(xt + xJ) and one bosonic band of width Wh = 

From analogy to the conventional t-J model the slave- 
boson representation with its fermionic description of 
low-lying orbital excitations is expected to describe well 
an orbitally disordered state far from critical instabilities 
towards orbital ordering. The existence of such an or- 
bital liquid state is strongly supported by the analysis of 
Ishihara, Yamanaka, and NagaosaQ who find orbital dis- 
order even within a slave-fermion Schwinger-boson rep- 
resentation which favors orbital ordering. Assuming the 
stability of the disordered state, we select the slave-boson 
representation, noting, however, that short-range orbital 
correlations will not be fully captured in this picture. 

The memory function, Eq. (||), is calculated via the 
force-force correlation function, replacing holon opera- 
tors by their mean-field value and keeping fluctuations 
around this value up to first order. Diagram (a) of 
Fig. [I] describes transitions between the two coherent 
Gutzwiller bands of the spinon mean-field Hamiltonian. 
These interband transitions are allowed due to the fact 
that the t term of Hamiltonian (|l|) is not diagonal in 
orbital quantum numbers. Using a free-fermion pic- 
ture, similar transitions were found by Shiba, Shiina, 
and TakahashitJ to contribute to the optical conductivity 
spectrum. The presence of electron correlations, how- 
ever, strongly suppresses the spectral weight associated 
with these processes and shifts the upper bound of the 
corresponding absorption spectrum to low energies. 

Diagrams (b) and (c) of Fig. |l| purely originate from 
correlations among electrons and are absent in a free 
theory. They describe transitions between a coherent 
Gutzwiller and a highly-incoherent band with corre- 
sponding contributions to the optical conductivity spec- 
trum that extend approximately up to bare fermionic 
bandwidth. We note that such broad incoherent ab- 
sorption was also found in a recent finite-temperature 
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(a) (b) (c) 

FIG. 1. Contributions to the force- force correlation func- 
tion. Solid, wavy, and dashed lines denote spinon, holon, and 
pseudospin Green's functions, respectively. 

diagonalization study performed on small clusters.El Di- 
agram (b) accounts for the composite nature of charge 
carriers in a strongly correlated system reflected by the 
spinon-holon convolution. The force operator employed 
in the calculation of this diagram yields terms involving 
two as well as three lattice sites. The two-site terms are 
already contained in the current operator and thus also 
appear in a direct calculation of the optical conductivity 
via the current-current correlation function. The three- 
site terms, however, are intrinsic to the force-force cor- 
relation function and partially account for vertex correc- 
tions. Diagram (c) finally describes scattering of carriers 
on pseudospin fluctuations. We note that the force term 
-F/ neglected before contains only terms involving three 
lattice sites of small amplitude, supporting our previous 
approximation to discard this term. 

The spinon and holon propagators of Fig. [I] are calcu- 
lated with the self-energy corrections shown in Fig. [|. 
The spinon self-energy diagrams (a) and (b) describe 
scattering on pseudospin and holon-density fluctuations, 
respectively. The quasiparticlc weight, which is propor- 
tional to x on the mean-field level, is further suppressed 
by these incoherent scattering processes. This small- 
ness of quasiparticle weight is consistent with the ex- 
perimental observation of only a small discontinuity in 
the spectral-weight functk)n-.at the Fermi energy as seen 
in photoemission spectra.H~El The bosonic self-energy di- 
agrams (c) and (d) describe scattering on fluctuations of 
the spinon bond-order parameter, creating a broad inco- 
herent peak in the bosonic spectral function. 

The expressions obtained from the diagrams in Figs. ^ 
and |^ contain several integrations over momentum space 
which we solve numerically using a Monte Carlo algo- 
rithm. The zero-frequency current-current correlation 
function xo of Eqs. (g) and (||) is obtained from the 
spinon mean-field Hamiltonian as xo — 2xxte 2 /a. We 
choose a hole-doping concentration of x = 17.5% and 
set J — QAt. The value of the spinon bond-order pa- 
rameter is numerically determined to be x = 0.25, the 
lattice constant is set to a — 3.9 A. The real part of the 
optical conductivity, Eq. (||) , consists of the Drude com- 
ponent ctd(w) = ttXo [l + 3fJ(rfM/dw)| a; ^o] 1 $( UJ ) and the 
regular part a reg (uj). The spectrum of the regular part is 
shown in Fig. ^| indicated by a solid line. We compare our 
theoreticaLresult to experimental data obtained by Oki- 
moto et alEfor 17.5% doped Lai^Sr^MnOa at T = 9 K 
represented by the dashed line in Fig. ||. The only free 
parameter which we use to fit the theoretical to the ex- 
perimental curve is the free fermionic bandwidth which 




(c) (d) 
FIG. 2. Self-energies of (a), (b) spinons and (c), (d) holons. 

we fix by setting 2-jj= 0.36 eV consistent with band struc- 
ture calculations.^ 

Good agreement between experiment and. theory is 
found for intermediate and high frequencies. Ej We stress 
that no additional fitting parameter is needed to obtain 
the correct absolute values of cr reg (w). The total spectral 
weight consisting of the Drude part ud{u) and the inco- 
herent part a reg (uj) is in agreement with experiment as 
can be seen from the effective charge-carrier concentra- 
tion defined as 

N e g = 2to °° I duj(<j D (uj) + <7 Icg (uj)); (6) 
ne Jo 

theory and experiment]^ both yield N c g — 6%. 

In the low-frequency region the theoretical curve devi- 
ates from the experimental one as a pseudogap opens in 
the spectrum. As a result our theory does not completely 
account for the weight transferred to the incoherent part 
of the spectrum as is reflected by the ratio of total spec- 
tral weight to Drude weight; here we obtain a value of 
2.8 which compares to « 5 found experimentally.!] 

The discrepancy between theory and experiment in the 
low- frequency part of the optical absorption spectrum in- 
dicates an additional scattering mechanism to be active 
at low energies which is not incorporated in our theory. 
We speculate this mechanism to stem from the closeness 
of the real system to orbital ordering underestimated in 
the above slave-boson treatment. Scattering on low- lying 
orbital collective modes induced by electronic superex- 
change and Jahn- Teller coupling as well as scattering on 
phonons will enhance the low-energy region of the spec- 
trum, thus filling the pseudogap. To make the latter 
point more explicit we calculate the conductivity includ- 
ing electron-phonon effects. Coupling of e g pseudospins 
to double-degenerate Jahn- Teller phonons and of charge 
to the lattice breathing mode are described by 

ffjT = -Aj T J2 (K ! + a *K + *04 - a *K) > ( 7 ) 

i 

H ch = -A ch ( b l + b i) (1 - rii) , (8) 

i 

respectively. For simplicity, phonons are considered to 
have dispersionless energy, which we set to loq = 0.05 eV. 
Corrections to the force-force correlation function and to 
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manganites. This transition, which due to the double- 
exchange mechanism is accompanied by a reduction of 
the mobility of charge carriers, suppresses the energy 
scale (oc xt) of orbital fluctuations and enhances the 
Jahn- Teller couplingia of these fluctuations to the lat- 
tice. Apparently, the picture of an orbital liquid which 
is quantum disordered by the metallic motion of holes is 
no longer valid above T c ; a disorder-order crossover in 
the orbital sector is in fact indirectly indicated by the 
temperature dependencapof local lattice distortions ex- 
perimentally observed. EZllia We speculate this crossover 
to control the metal-insulator transition by further re- 
ducing the mobility of charge carriers, thus playing an 
essential role in enabling the formation of small lattice 
polarons and their localization above T c . 

We wish to thank P. Fulde and P. Horsch for valuable 
discussions. 



FIG. 3. Incoherent part of the optical conductivity <j tbs (uj) 
as a function of photon frequency lu. The solid line repre- 
sents the theoretical curve for x = 17.5%, which is fitted ta 
experimental data (dashed line) obtained by Okimoto et alB 
for Lai-ajSr^MnOa at T — 9 K. The dot-dashed line corre- 
sponds to the theory including lattice effects. 



the fermionic self-energies are evaluated within a weak- 
coupling scheme, assuming the dimensionless coupling 
constant C = ^ 2 N(ep)/uJo, where N(ep) is the total e g - 
density of states at the Fermi level, to be small below T c . 
The result, which for Qt = Cch = 0.3 is shown by the 
dot-dashed line in Fig. ||, suggests that (possibly strong) 
lattice effects are present even in the metallic ferromag- 
netic phase. 

We further calculate the constant of T-linear specific 
heat 7 from the spinon mean-field Hamiltonian. For 
x = 17.5% we find 7 = 7.2 mJ/mol K 2 as compared 
to experimental values 5-6 mJ/mol K 2 (Ref. ^) and 
3.3 mJ/mol K 2 (Ref. ^). Experimentally 7 was observed 
to be nearly independent of x. In our theory we find 
7 oc tN°(ep)/(xt + xJ)', since the bare density of states 
iV°(e) exhibits a pseudogap centered around the chem- 
ical potential at half-filling, 7 is rather insensitive to 
changes in x for moderate hole-doping concentrations; 
for x — 30% we in fact find 7 = 7.1 mJ/mol K 2 . 

In summary, we have calculated the optical absorp- 
tion spectrum of ferromagnetic manganites, emphasiz- 
ing the role of low-energy orbital fluctuations leading to 
the strongly correlated nature of the metallic state. The 
theory explains the large incoherent component of the 
optical conductivity accompanied by a strong suppres- 
sion of Drude weight observed experimentally and is also 
consistent with measurements of the specific heat. The 
fact that the anomalous transport properties in the fer- 
romagnetic, phase are described well supports the orbital- 
liquid ideaB and suggests that the orbital degrees of free- 
dom coupled to the lattice are of relevance in the metal- 
insulator crossover driven by the magnetic transition in 
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